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In superfiuid 3 He turbulence is carried predominantly by the superfiuid component. To explore 
the statistical properties of this quantum turbulence and its differences from the classical counter- 
part we adopt the time-honored approach of shell models. Using this approach we provide numerical 
simulations of a Sabra-shell model that allows us to uncover the nature of the energy spectrum in the 
relevant hydrodynamic regimes. These results are in qualitative agreement with analytical expres- 
sions for the superfiuid turbulent energy spectra that were found using a differential approximation 
for the energy flux. 
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Introduction: Supcrfluids are known to display many 
spectacular properties, such as the fountain effect and 
rolling films, which arc unique to quantum hydrodynam- 
ics. Ever since the discovery of superfiuid it was con- 
jectured that their hydrodynamic can be described as 
two inter-penetrating fluids with two different velocities 
v s (superfiuid component) and v n (normal component )i. 
The relative density of each component, (p s and p„ re- 
spectively), depends on the temperature; this so-called 
"two-fluid" model was later developed very usefully by 
LandavJ.. In this model the normal component is de- 
scribed by the Navier-Stokes equations and the superfiuid 
part is modeled by the Euler equation. Accordingly the 
inviscid super-component is characterized formally by an 
infinite Reynolds number: it is very easy for superflu- 
ids to become turbulent. Quantum mechanics plays an 
imprtant role: unlike classical fluids where the circula- 
tion around the vortices is a statistical and dynamical 
variable, the properties of quantum vortices are tightly 
constrained by the laws of quantum mechanics. The cir- 
culation around a quantum vortex is quantized to integer 
values of a fundamental unit called the quantum of circu- 
lation k = 2-Kh/m, where m is the mass of a superfiuid 
particle. Because of this one can identify a new typical 
length-scale in quantum turbulence, refereed to herewith 
as £, which corresponds to the mean distance between 
quantum vortices. 

At finite temperatures T ^ the normal and super- 
fluid components are coupled to each other via thermal 
excitations and the resulting dissipative force is called 
the mutual friction. This description is in qualitative 
variance with ordinary fluids and accordingly the ques- 
tion of how different classical and quantum turbulence 
are has emerged as a hot topic in the fluid mechan- 
ics community^. An important step in establishing a 
theory of quantum turbulence consists in developing a 
framework that describes the energy spectra of both the 
normal E n (k) and the superfiuid E s (k) components. The 
work to achieve this goal is still developing and there is a 
possibility that progress here might also shed some light 
on some long standing classical turbulence problems. 

Because of the high kinematic viscosity of 3 He (com- 
parable to that of olive oil)*, the normal component can 
usually be considered to be at rest with respect to the 
superfiuid component with v n <C v s . In such circum- 



stances we propose that the overall behavior of quantum 
turbulence is largely dominated by the superfiuid veloc- 
ity v = v s component only. Therefore, the description of 
the fluid motion can be reduced to the coarse-grained hy- 
drodynamic equation for the superfiuid velocity v(r,t), 
considered on length scales of the order of I which are 
larger than intervortex distance £: 
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Here p is the chemical potential, f(r,t) is a random force 
which is introduced in order to mimic the energy injec- 
tion that sustains turbulence; this force acts at the outer 
scale Iq ^> I ^> i. The dissipative term D represents 
the effective mutual friction which is the last remnant of 
the underlying but vanishingly small normal component. 
Throughout the paper it should be kept in mind that the 
largest relevant wavevector is fc max = in other words, 
we restrict ourseleves to the pure hydrodynamic regime 
of superfiuid turbulence. In this case, the mutual friction 
can be approximated* as D w a u) X (u> x v) /|w|. Here a 
is the temperature dependent dimcnsionlcss mutual fric- 
tion parameter and u) = V x v is the superfiuid vorticity. 
This expression can be further simplified by averaging 
out the vectorial structure and recognizing that vorticity 
in hydrodynamic turbulence is dominated by the small- 
est possible eddies with k ~ l/l. Because of their small 
size, these eddies also have very small turnover times. 
On the other hand, the superfiuid velocity v is domi- 
nated by the largest eddies who, in turn, have a very 
long turnover time. This means that u) and v can be 
considered as almost completely uncorrelated. In this 
case, we can replace uj by its mean square value and the 
dissipation term can be approximated by 



D w Tv, 



T = acj T , c4 = (|u;| 2 )«2 / 

Jko 



(lb) 

k 2 E(k) dk , (2) 



where E(k) is the so-called 1-dimensional energy spec- 
trum, normalized such that the total energy density per 
unite mass is E = J E(k)dk. 

Our goal below is to analyze the stationary energy 
spectra of superfiuid turbulence, E(k), at scales larger 
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than the intervortex distance. This allows us to sim- 
plify further the problem by considering T in Eq. (jlb|) 
as a time independent external parameter. Then, after 
finding E(k) for a prescribed T we can get the a corre- 
sponding to that T using Eq. @. We will see that even 
after all these simplification the spectra found below are 
highly non-trivial, justifying this study on its own right. 

In Sec. I we examine Eqs. (TTJ) using a shell model of 
the type that was widely used in studies of classical tur- 
bulence. Of course, the shell model is an uncontrolled 
approximation of Eqs. (H}. Nevertheless our experience 
with classical turbulence is that they provide strong in- 
dications to the correct spectra associated with the full 
equations. This Section presents numerical simulations 
that resolve the energy spectra E SM (k) for different values 
of the mutual friction extending over an unprecedented 
range of wave vectors; see the left panel of Fig. [T] We 
demonstrate that energy spectra depend crucially on the 
value of dimensionlcss mutual-friction frequency 
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which is normalized by the energy influx eg at the outer 
wave vector ko ■ There is a critical value 7 cr sa 1 at which 
E(k) = E CT oc fc~ 3 . For 7 < 7 cr one has supercritical 
spectra E sp (k) > E CT (at the same values of £0 and ko) 
with Kolomogorov- 41 (K41) tail E sp (k) ~ e 2 / 3 fc" 5 / 3 and 
smaller energy flux e < sq. For 7 > -f cr one has subcriti- 
cal spectra E s ^(k) < E cv that terminate at final k. 

To increase our trust in this complex behavior we re- 
analyze (in Sec. II) the same basic Eqs. (TTJ) analytically, 
using approximate differential closures for the energy flux 
over scales. These studies culminate with analytical so- 
lutions Eqs. (fTUj) and (fl"3l) for the super- and subcritical 
energy spectra respectively. The closure techniques are 
not exact cither, but we will see that they agree very 
well with the results of the shell model. The analytical 
spectra for E sp (k) and E s b(k) are shown in the middle 
panel of Fig. [1] where the stars correspond to the numeri- 
cal data from Sabra-shell model with similar values of the 
ratio 7/7cr- One sees very good qualitative agreement be- 
tween numerical and analytical results, despite the fact 
that both approaches are (uncontrolled) approximations 
of a very different nature. We propose that our results 
for the energy spectra can be considered as valid predic- 
tions on a semi-quantitative level. A detailed analysis of 
our results is presented in the final Sec. III. 

I. Superfluid energy spectra in shell models 

Here we investigate the superfluid energy spectra nu- 
merically in the shell-model approximation. Compet- 
ing methods, such as vortex filament simulations or the 
solution of the non-linear Schrddingcr equation are ex- 
tremely demanding computationally, having therefore a 
very short span of scales. Shell models offer an at- 
tractive alternative, as recently suggested by Wacks and 
Barenghi^ who proposed a two-fluid shell model for su- 
perfluid turbulence using a GOY shell model^ with ad- 



ditional coupling terms. While only 18 shells were used 
in Ref. Q, we have been able to expand the range of 
scales quite significantly by increasing the number of 
shells to 36. This improvement allows us to probe an 
inertial interval over scales which extends to more than 
five decades, thereby making it possible to resolve sub- 
regions with different scaling behavior. Also, instead of 
using a GOY model, which suffers from the artifact of 
"period-three" oscillations in the energy spectrum, we 
used the Sabra shell model^ in which this non-physical 
phenomenon is removed and the scaling behavior become 
more transparent. These improvements allow us to get 
essentially new information about the statistical behavior 
of superfluid turbulence. 

A. "General" shell models of hydrodynamic turbulence 

Shell models of turbulence£&i£~— are simplified cari- 
catures of the equations of fluid mechanics (jla[) in wave- 
vector representation: they mimic the statistical behav- 
ior of fc-Fourier components of the turbulent velocity field 
u(k,t) in the entire shell of wave- vectors k rn < fc m +i by 
only one complex shell velocity u m . The integer index m 
is referred to as the shell index and the shell wave num- 
bers are chosen as a geometric progression k m = fcoA m , 
with A > 1 being the shell spacing parameter; often one 
chooses A = 2. The equation of motion of the fluid [e.g. 
Fourier transform of Eq. (|laj) ] is reduced to a dynamical 
equation of motion for u rn (t): 



+ r 



u m = NL m {wj} + /„ 



(4a) 



The model contains a time-independent mutual friction 
parameter T and random force terms f m whose contribu- 
tion as the energy influx at large scales is usually limited 
to the first few shells only. The nonlinear term NL m {uj}, 
which originates from the (v-'V)v in Eq. (JTa)), should be i) 
proportional to k m (because of V), ii) quadratic in the set 
of velocities {uj} and iii) should conserve the total kinetic 
energy. In order to implement Richardson-Kolmogorov 
idea of a step-by-step cascade energy transfer over scales, 
the set of the shell indexes j in NL m {uj} should be local- 
ized around m. As as rule, one takes j = to, mil, mi2. 

B. Sabra-shell model of turbulence 

It is known that the original GOY -shell modeU& suffers 
from non-physical "period-three" oscillations of the en- 
ergy spectrum that originate from artificial phase corre- 
lations in subsequent triads of the shell velocities. These 
oscillations were eliminated in the Sabra-shell model — by 
a proper choice of the nonlinear term: 

NL m {u-,} = i(a k m+1 u m+2 u* m+1 (4b) 
+ bk m u m+ iu* n _ 1 - cfc m _iu m _iu m _ 2 ) , 

where * stands for complex conjugation and where the 
coefficients a, b and c are real. Conservation of energy in 
the forceless inviscid limit requires that a + b + c = 0. It 
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FIG. 1: Left/Middle panel: Numerical results for the energy spectrum (left) and the energy flux (middle) in the Sabra-shell 
model for different values of 7. Right panel: Analytical trial functions E sp (k), Eg . (|10|1 . and E s b(k), Eg. (|13|l . in the differential 
energy- flux model. The stars show a comparison with data from the Sabra-shell model at similar values of the ratio 7/701 ■ The 
bottleneck effect E Bp (k) ~ k 2 at the highest wavevectors is removed from the Sabra-shell model by proper choice of parameters 
to simplify the simulations. 



is knowr>£~— that the statistical properties of the veloci- 
ties in shell models mimic those of the three-dimensional 
hydrodynamic turbulence for a = 1, b = c = — 0.5. In 
the forceless inviscid limit the model has two quadratic 
invariants: the energy E and H, often associated with 
helicity: 



E E 



(5a) 



The ID turbulent energy spectrum E(k) can be con- 
nected to E m by the relation E(k m ) = (E m )/k m . The 
energy flux over scales e m in the Sabra model is: 



e(k m ) = 2lm[ak m+ iS m+ i - ck m S m ] 



(5b) 



where S m 
tion function and Im[ 



is the third-order correla- 
denotes imaginary part of [. . . ]. 



C. Numerical realization of the Sabra model 

The set of N = 36 coupled Eqs. (j4|) was solved us- 
ing 4th order ETD-Runge-Kutta mcthocU 3 - for stiff sys- 
tems. In order to mimic the energy sink at very high 
wavenumbers (e.g. due to Kelvin waves), a hyperviscos- 
ity term of the form vk 4 \u m \ 2 with v = 10~ 32 was added 
to the equation. Energy pumping into the system was 
carried out by the "helicity-free" forcing f„JL. We used 
random forcing (5-correlated in time ), with a finite am- 
plitude and a random phase in the first two shells, such 
that th e forcing amplitudes are related to each other as 
f-2 = sf— (c/o)/i and f\ = 0.01. Such a forcing mini- 
mizes the input of helicity into the system. The mutual 
friction coefficient 7 was chosen in the range < 7 < 1.65 
to expose the different solutions exhibited by the system. 
The use of hyperviscosity allowed to maximize the extent 
of the inertial interval. However the sharp onset of the 
viscous term lead to artificial oscillations in the energy 
spectra in a few shells around the dissipative cutoff (see 
Fig HI left). These oscillations are an artifact of the nu- 
merical procedure and not the result of the underlying 
dynamics. 



Results of the numerical simulations of Sabra Eqs. (j4|), 
shown in Figs. [1] and [3J will be discussed below in Sec. III. 

II. A differential model for the energy spectra 

A. Energy balance equation with differential closure 

The energy spectrum of pure hydrodynamic turbulence 
can also be studied analytically by investigating the en- 
ergy balance equation 



deQe) 
dk 



TE(k) =0, 



(6) 



which is presented here in the stationary case 
(i.e. dE(k)/dt = 0). The energy density E(k) is related 
to the total energy E by: E = J E(k)dk. Equation © 
directly follows from Eqs. ([1} with some expression for 
the energy flux over scales, s(k), in terms of the third 
order correlation function of Uf, exactly as in the case of 
classical turbulence. In order to proceed further one can 
borrow a closure procedure from classical turbulence that 
expresses e(k) in terms of the energy spectrum, E(k). 
Even though this step is widely used, it is worth remem- 
bering that it is uncontrolled. The simplest algebraic 
closure relation suggested by KovasznayJ^ 



e{k) ~ ^ 2 E 3 ^ 2 ; 



(7a) 



just follows from K41 dimensional reasoning. The pref- 
actor I is chosen to simplify the appearance of some of 
the equations below. Notice that this simple algebraic 
approximation ()7a[) fails to include an important effect: 
the bottleneck energy accumulation near the inter- vortex 
scale I. We therefore use the Leigh differential closure^ 
for the energy flux in the form suggested by Nazarenko-^: 



1 



s(k) = - 8 Vk^E- 



E(k) 



k 2 



(7b) 



Dimensionally, this closure relation coincides with 
Eq. d7aj), but replaces E(k)/k 3 with d[E(k)/k 2 ]/dk. This 
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allows one to account for the existence of thermodynamic 
equilibrium when e(k) = and E(k) oc fc 2 . The prefac- 
tor — | is chosen to get a frictionless (i.e. with T = 0) 
K41 solution, 

E K41 (e\k) = C K e 2 / 3 k- 5 /\ C K = (^f'* « 1.68, (8a) 

with C K numerically close to its experimental value. 

B. Simplified solution for the energy spectra 

In this Section we overview some simplified solu- 
tions for the energy spectra that will serve as build- 
ing blocks for the general solution. An exact solu- 
tion of Eqs. dp with the differential closure (|7b[) was 
found in Ref. in the case V — 0. We rewrite it as 
E{e\k) = E KA1 (e\k)T cq (k) with 



T cq (fc) = [l+ (fc/fc e 



ll/2-i 2/3 



(8b) 



One can sec that E(e\k) has a thermodynamic tail 
E{k) oc k 2 for k 3> fc eq , as required. A preliminary anal- 
ysis into how to account for the mutual friction term in 
the balance equation Eq. (j6]) was done in Ref. [|| by us- 
ing the algebraical closure (|7ap and three regimes were 
identified: 

• "Critical" scale-invariant solution 



E CI (k) 



16 £2 
25 F 
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k 



fc :S 



(9a) 



At the special value 7 = 7„- = ^-q— ~ 1.46 it is realized 
for any k. One sees that 7 cr is of the order of unity 
and that the corresponding dimensional T cr is about the 
turnover frequency of the outer-scale eddies. 

• For 7 > 7 cr solutions become "sub-critical" £" s b (k) < 
E cv (k) and vanish at some finite k = fc cr : 

E sh (k)~E cr (k)T„(k), T_(fc) = (l- ^J) 2 . (9b) 

• When 7 < 7 cr the energy spectra become "super- 
critical" with E s b(k) > E CT {k). For k -C k CT , these spec- 
tra are close to the critical one and for k S> fc cr they 
develop a K41 tail parametrized by < Sq- 



E sp (k) ~ £ K41 ( £oo |fc):r + (fc), T+(k) = 1 



1. 2 / 3 

"^cr 
k 2 ^ 



(9c) 



When 7 — > 7 cr , — > and fc cr — > 00: 

£00 =£o(7-7cr) 2 , fccr = fco|7^7crr 2/3 • (9d) 

These simple solutions allow us to construct physically 
motivated approximate analytical solutions to Eqs. ([5J 
I7b[) . The accuracy of these solutions can be improved if 
necessary by introducing a correcting polynomial expan- 
sion to our trial energy spectra. However, for the purpose 
of the present discussion it is sufficient to present only 



the most simple versions of the energy spectra which are 
determined by physical considerations only. 

Notice that with the differential closure (|7bp the energy 
balance Eq. ^ has a scale invariant critical solution (|9"aj) 
at precisely the same value of 7 cr . Moreover, for 7 
7 cr there are sub- and super-critical solutions, E s \,(k) < 
E cr (k) and E sp (k) > E CI (k) that we will discuss below. 

C. Supercritical energy spectra 

Detailed analysis of the differential energy balance 
Eqs, ((HI I7bp leads us to the following form of approxi- 
mate super-critical energy spectrum for T < T CT : 



E sp (k) 
T sp (fc) 



E Kil (s\k) { ^T 2 (k) + T 2 q (k)-1 - T sp (k) } , 

"ii\!/3 r / k 

15 ' 



(ek 2 )V3 Vfc, 



7/6 



(10) 



This solution is plotted in Fig. [TJ middle panel. The 
solution (fTO)) is constructed from the building blocks 
-E'K4i(£|fc) i T eq (k) and 7+(fc) discussed above and given 
explicitly by Eqs. (|5a|) . (|8b| and (JScJ) respectively. This 
solution (|10p has a simple physical interpretation, which 
will be clarified by spelling out the regions of k which 
possess different scaling behavior; these regions are sep- 
arated by fc cr and fc eq : 

1, T+ « (fc cr /fc) 4 / 3 » 1, T sp « 1: 



• For k < fccr, T cq 



.2/3L4/3 

E sp (k)^C K £ fc 3 , e(k) 



i.e. solution has the form close to critical. 

• For fccr < fc < fc oq : T+(fc) a T cq (fc) » 1, T sp (fc) < 
1. Thus in this region one observes the K41 spectrum 
E K41 (e\k) with the energy flux e that is related to fc cr 
and r by the following global constraint 



£0 



# sp (fc)dfc 



(12a) 



Two other relations between these parameters are found 
from the boundary conditions for E sp (k) and e{k) at fc = 
fco- For example from e(fco) = £0 one finds: 



£f) 



15 / fccr y 
11 V fc n J 



(12b) 



Notice that square root in Eq. (|10p is constructed such 

that for fc > fc cr + T 2 q ~~ 1 w T oq! thus ensuring 

that E sp (k) is close to the exact solution of Eqs. (0 I7bp 
at r = 0, that predict the bottleneck energy accumula- 
tion for fc > fc cq . This wavevector should be fixed by 
the boundary conditions at smallest possible scale in the 
problem, inter-vortex distance I. This condition is out- 
side of the scope of the present description and currently 
we should consider fc oq as an external parameter. 

Finally, the T sp (fc) term was obtained by performing 
an asymptotic analysis of the correction to the solution 
(flU when f^O and fc > fc oq . 
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FIG. 2: Left/Middle panels: Mismatch functions Xsp (resp. Xsb) for the supercritical trial solution (resp. subcritical) defined 
in the text. Right panel: Evolution of the position k CI as a function of the strength of the dissipation parameter. The vertical 
dashed line corresponds to T = F CT and the horizontal dashed line represents fc cr = feo- 



The accuracy of the trial solution E sp (k) can be quan- 
tified by introducing a mismatch function Xsp(k) = 

100 x [ de ^ /TE sp (k) — 1] which measures the percent- 
age of error of E sp (k) when injected into Eq. ((HJ). As we 
can see in Fig. [5] (left panel) , the disagreement never ex- 
ceeds 15% except in the vicinity of fc oq which has been 
left as an outside parameter here. Having in mind the ap- 
proximate character of the differential approach we con- 
clude that our analytical form (|10p for E sp (k) is a good 
approximation to the exact solution. Recall also, that 
if needed, our solution can be improved by supplement- 
ing it with appropriate polynomial corrections in the way 
demonstrated in Ref. (l7j |. 

B. Subcritical energy spectra 

As 7 — > 7 cr , the supercritical spectrum smoothly con- 
nects with the critical solution E cr (k) and indeed fc cr — > 
+oo as can be seen in FigfTJ middle panel. Eventually, 
we reach a point where 7 > 7 cr and the dissipation is so 
strong that the differential approximation predicts the 
resurgence of a finite wavevector fc cr such that the en- 
ergy spectrum E s i>(k cr ) = as well as the energy flux 
s(fcsb) = vanish. In this case, we propose the following 
form of approximate sub-critical energy spectrum 



E sh {k) = E CI (k)T_(k)T sh (k) , T sh (k) = 



1 - 



k 3 / 4 



k 



3/4 



(13) 



This solution is plotted in Fig. [TJ middle panel. In addi- 
tion to the building blocks already discussed above, this 
expression contains an extra correction T s b(fc), that fixes 
its asymptotical behavior close to the terminal wavevec- 
tor fc cr . Conservation of energy, in the form of the global 
constraint 



fccr 



E B (k)dk = £ , 



(14) 



selects the value of fc cr , whose evolution as a function of T 
can be seen in Fig(2J middle panel. When T — > T CT the po- 
sition of the terminal wavevector is pushed to k ct — > +00 
which is consistent with a smooth connection with the 



critical solution E cr (k). On the other hand, when the 
dissipation is very strong T ^> r„, the terminal wavevec- 
tor rapidly decreases and approaches k €I — > fco asymptot- 
ically as it should. 

The accuracy of the trial solution E s b(k) can be as- 
sessed by defining a mismatch function Xsb in complete 
analogy with that defined in the previous paragraph. We 
can see in FigJ5] (middle panel) that the error always stays 
below 12%. This observation leads us to the conclusion 
that E s b(k) is a good analytical approximation to the 
exact solution which is more than satisfactory for our 
purposes. 

III. Discussion and conclusive remarks 

A. Numerical vs. analytical energy spectra 

The main results of our numerical simulations in the 
framework of the Sabra-shell model are presented in 
Fig.Q] In order to continue our discussion, notice that in 
the left panel one sees a horizontal line for 7 = in the 
plots of k 5 / 3 E(k) vs. k. Clearly, a K41 energy spectrum 

E(k) tx £Q^ 3 fc -5 / 3 and fc-independent energy flux e(k) 
which is equal to the energy influx £0 (shown in right 
panel) is expected in this case. Intermittency corrections 
are practically invisible. 

For small values of mutual friction parameter 7, we see 
that the energy flux first decays in the region of small k 
before approaching £00 whose value is analytically esti- 
mated by Eq. (|9d[) . This behavior can be understood 
quite straightforwardly. Indeed, while the mutual damp- 
ing frequency T is fc-independent the K41 eddy-lifetime 
frequency increases as w(fc) ~ e 1 / 3 ^ 2 / 3 . Therefore it is 
always the mutual damping that must dominate the en- 
ergy flux over scales in the region of small k. Such a 
damping naturally results in a decay of £(fc). On the 
other hand, in the region of large k it is the characteris- 
tic frequency of energy transfer over scales, uj(k), which 
dominates over T and therefore the mutual friction can 
be neglected. In this region e(k) — > £oo and the K41 scal- 
ing E(k) ' must be recovered as is indeed confirmed 
both by the numerical and analytical plots of E(k) in the 
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left and middle panels. 

Moreover, physical intuition [supported by Eq. (|9d|) ] 
tells us that for some (critical) value of 7 the asymptotical 
value limfc-j.oo = should vanish. In this critical regime 
the mutual damping frequency T is balanced by the cas- 
cade frequency ui(k) which can be estimated as ky / kE(k). 
This immediately suggests that E CI (k) ~ T 2 /k 3 which 
is corresponding precisely to our numerical observations 
(dashed red line in left panel) as well as our analytical 
predictions, Eq. (|9"aj) . 

Finally, we conclude by mentioning that the numeri- 
cal value of fc cr is expected to be very sensitive to 7 in 
the vicinity of 7 cr as hinted at by Eq. (|9d|) . Indeed, the 
complete evolution of k cr in both super- and subcritical 
regimes is presented in Fig[2] (right panel) . 

B. Mutual friction parameters 7 vs. a 

The transition from subcritical to supercritical energy 
spectra is controlled by the dimcnsionles parameter 7 in- 
stead of the more experimentally accessible mutual fric- 
tion parameter a. Both can be related to each other by 
Eqs. © and Notice however that unlike -f CI ~ 1 

which is universal, the critical value a cr does depend on 
experimental parameters such as the energy injection rate 
at the outer scale I and on the extent of the "classical" 
inertial interval l/t. It is nevertheless interesting to esti- 
mate a cr on the basis of typical experimental parameters. 
Substituting the critical solution ((9a)) in Eq. ([2} one gets 



" 4 v /2hi(fc ^) ^Ry) ' 

With k e ~ 10 3 [3 this gives a cr ~ 0.34. For a < a cr 
one can neglect the effect of the mutual friction and di- 
rectly substitute the K41 spectrum (j8aj) into Eq. (|2|). 
Under the assumption that this spectrum is valid up 



to k ~ l/£, one finds a remarkably simple relation for 
small a: 

7 -2.25 a. (16) 

One should be mindful that due to the possible bottle- 
neck energy accumulation^, which depends on the ra- 
tio of £ to vortex-core radius, the numerical coefficient 
in Eq. (|16[) can increase. Notwithstanding these reser- 
vations, we believe that this equation should be useful, 
even if only as a rough estimate. 

C. The road ahead 

To finish the discussion we should underline again the 
approximations that will have to be relaxed in future 
work. 

• The first one is the neglect of the turbulent motion 
of the normal component. We can study this approxi- 
mation numerically, using a properly generalized^ Sabra- 
shell model^, and analytically, within properly formu- 
lated differential approximation for the normal and su- 
pcrfhhd energy fluxes. 

• Secondly, using the two-fluid results we will gener- 
alize the model of Ref of gradual eddy-wave crossover 
to the case of non-zero temperatures, to study the ef- 
fect of temperature suppression of the bottleneck energy 
accumulation near intervortex scales. 

• Lastly we did not consider stability; the analysis of 
the stability of one- and two-fluid energy spectra against 
thermal perturbations should be accomplished, having in 
mind the possible self-sustained oscillations around, or in 
the vicinity of stationary solutions. 
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